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ABSTRACT 


The MSFC bearing seal material tester (BSMT) can be used 
to evaluate the SSME high pressure oxygen turbopump (HPOTP) 
bearing performance. The four HPOTP bearings have both an 
imposed radial and axial load. These radial and axial loads 
are caused by the HPOTP' s shaft, main impeller, preburner 
impeller, turbine and by the LOX coolant flow through the 
bearings respectively. These loads coupled with bearing 
geometry and operating speed can define bearing contact 
angle, contact Hertz stress and heat generation rates. The 
BSMT has the capability of operating at HPOTP shaft speeds, 
provide proper coolant flowrates but presently, can only 
apply an axial load. Due to the inability to operate the 
bearings in the BSMT with an applied radial load, it is 
important to develop an equivalency between the applied 
axial load and the actual HPOTP loadings. 

In this study, the objective was to use the 
SHABERTH/SINDA ( shaf t-bea r ing- the rmal ) computer code to 
simulate the BSMT bearing-shaft geometry and thermal-fluid 
operating conditions. This study was performed at two shaft 
speeds using two coolants, LN2 and LOX. A simulation of the 
HPOTP was also generated by SRS/System Division using 
current operating conditions from the SSME HPOTP. Then, a 
comparison of the bearing contact stresses and heat 
generation rates of these two simulations was attempted to 
establish the equivalence between the BSMT axial load and 
the HPOTP loads. 
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INTRODUCTION 


In the Space Shuttle Main Engine (SSME) High Pressure 
Oxygen Turbopump (HPOTP), four ball bearings support a 
turbopump shaft, a main impeller, preburner impeller and 
turbine. Throughout the flight history of the SSME, these 
bearings have been subject to various degrees of damaginq 
wear. Two possible causes for this wear are insufficient 
lubrication resulting in frictional heat generation and 
large contact (Hertz) stresses between the balls and the 
inner and outer races due to loading and bearing geometry 
variations. Even though these causes will be addressed in 
this study, numerous scenario's based on test data can be 
formulated to address the HPOTP bearing wear problem. The 
mam source of test data is from instrumentation 
measurements of the HPOTP. However, due to the expense of 
this process, viable alternatives to predict bearing 
behavior must be established. One alternative is the use of 
the NASA-Marshall Space Flight Center (MSFC) Bearing Seal 
Material Tester ( BSMT ) . Another relatively inexpensive 
alternative is to develop a computer model to simulate the 
bearing environment. A general program called SHABERTH 
(Shaft-Bearing-Thermal) developed originally by SKF 
Industries and later greatly modified by SRS 
Technologies/System Division exists and will be used to 
attempt this simulation. In addition to SHABERTH which 
analyzes the bearings and shaft, a code named SINDA (System 
Improved Numerical Differencing Analyzer) will be coupled to 
SHABERTH to perform the temperature calculations. Thus, 
this code will be referred to as SHABERTH/SINDA. 

The major unknown in this study of bearing behavior is 
loading. From experimental studies on the HPOTP, Figure 1 
shows the best estimate of the loads applied to the shaft 
the preburner impeller, main impeller and turbine 
that the bearings support. In addition to these radially 

loads, there also exists axially applied loads due 
to the pressure«area (PA loads) of the liquid oxygen (LOX) 
coolant that flows through the bearings. These PA loads 
are of particular importance when the turbopump throttles 
its speed. 

Figure 2 shows a schematic of the bearing-shaft 
arrangement and the flow paths through the BSMT. To 
reproduce HPOTP conditions at this time is not possible 
since the tester has a different flow path than the HPOTP, 
the working fluid in the tester is LN2 (liquid nitrogen)' 
not LOX, and most importantly, there can be only an applied 
axial load in the tester to simulate PA loading and 
preloading. Thus, presently, no radially load can be 
applied to simulate the radial HPOTP loads. 
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TURBINE 



The purpose of this study is to attempt to use 
SHABERTH/SINDA programs to model the BSMT . This model will 
only have applied axial loads on the shaft and will be 
used in conjunction with a model of the HPOTP that was 
conducted by Spectra Research Systems ( SRS ) to compare heat 
generation rates and Hertezian stresses. Hopefully, this 
study will establish which applied axial loads for the BSMT 
model corresponds to the combined radial and axial loads for 
the SRS HPOTP model. From the comparison of heat generation 
rates and contact stresses, a so-called "equivalent" load 
can be stated for the BSMT based on HPOTP loading cases. 

Note that several important parameters as coolant flow rate, 
bearing geometry changes, coefficient of friction, coolant 
inlet temperature and pressure drop will be held fixed in 
this study. This was done to limit the problem's scope not 
to infer the insignificance of these parameter's affect on 
bearing behavior. in this study, only shaft speed will be 
varied along with type of coolant used ( LOX vs. LN2 ) . 

Recall, LOX is the coolant of the HPOTP, however LN2 is the 
current working fluid for the BSMT. The BSMT is currently 
undergoing redesign changes to eventually use LOX as the 
working fluid again. So, equivalent loads will be 
established using both fluids for the BSMT to simulate HPOTP 
loading . 


OBJECTIVES 

As previously stated, the purpose of this project is to 
simulate the BSMT conditions using the SHABERTH/SINDA 
computer code. Using this model of the tester and a 
turbopump simulation using SHABERTH/SINDA performed by SRS, 
a comparison of the heat generation rates and Hertz stresses 
will be made to attempt to correlate the axial load applied 
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Figure 2: The current BSMT and its flow paths 
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in the tester model to the axial-radial load combination 
that exists in the turbopump simulation. The objectives of 
this project were. 

1. To develop the input data necessary for 
modelling the BSMT using LN2 and LOX and 
perform a parametric study. 

2. To obtain SHABERTH/SINDA models of the 
turbopump from SRS/System Division. 

3. To compare for two different shaft speeds for 
both LN2 and LOX, the heat generation rates 
and contact Hertz stresses of two models to 
correlate the loadings applied to the tester 
simulation to those applied in the turbopump 
simulation . 


SHABERTH/SINDA Computer Models 

The SHABERTH program is structured in four sections: 
thermal, bearing dimensional equilibrium, shaft-bearing 
system load equilibrium and bearing rolling element and cage 
load equilibrium. A detailed account of these sections, 
bearing equations that are used, flowcharts of program 
structure, and sample input and output are described in 
reference (3). The bearing theory used in this problem is 
based on reference (1) by Harris. When SHABERTH was modified 
for the HPOTP by SRS , it was decided not to use the SHABERTH 
thermal model but to replace it with SINDA. SHABERTH uses 
an assumed set of temperatures given by a user then 
calculates all the bearing forces and moments, Hertz 
stresses, bearing geometry changes and heat generation 
rates. SINDA uses the calculated heat generation rates 
from SHABERTH to compute a temperature distribution. A 
UNIVAC computer runstream which controls the program flow 
replaces the assumed temperatures with the newly calculated 
SINDA temperatures. These temperatures that are being 
compared are of the shaft, inner ring, inner race, ball, 
outer race, outer ring, housing, bulk fluid temperature 
respectively. This iteration process between gHABERTH and 
SINDA continues until thermal convergence to 2 F occurs, or 
thermal runaway to 1000 F diverges the solution or when 15 
iterations occur usually related to an oscillating solution. 
Maximum runtime or maximum number of pages usually is 
associated with a divergence or oscillating solution. A 
good indicator of this type solution is when SINDA cannot 
reach an energy balance. For these cases of divergence, the 
SHABERTH/SINDA simulation will terminate. For convergence, 
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of the BSMT simulation 




the SHABERTH/SINDA simulation usually iterates 4 to 7 times 
depending on values of the initial temperatures assumed in 
SHABERTH. The computer run time of a converged solution is 
from 45 minutes to 1 hour. 

The input data to the SHABERTH model in general is 
discussed in reference (3). The Appendices to this 
reference are particularly helpful since it shows the 
formatting and structure of the input information and a 
listing of typical output. The input data that SRS added 
for their modifications to SHABERTH is described in (14). 
Much effort was expended to learn and verify analytically 
many of the inputs to SHABERTH and the source data in 
SINDA. However, some of inputs are based on experimental 
tester data. For instance, shaft dimensions and bearing 
locations, shown schematically in Figure 3, were found from 
the BSMT drawings. Fluid properties used for LN2 were 
found by interpolating at 480 psia, the tester pressure, 
using reference (5). In the same manner, fluid properties 
for LOX were found using reference (6). Cage load and 
viscous heat generation inputs were extensively calculated 
by myself based on J.C. Cody's notes from SRS Technologies. 
These calculations are based on the theory in reference 
(12). Cage heat generation rates based on the cage loads 
are found in a table in reference (14) as a function of 
coefficient of friction. 

In the Appendices of this report, a representative 
listing of SHABERTH input and references to the lines of 
SINDA code that are to be changed by the user are given for 
both LN2 and LOX. When shaft speed was varied, the inputs 
that must be varied were viscous heat generation rates for 
bearings 3 & 4 (VQBRGl, VQBRG2 ) , shaft speed (SHAFTS), cage 
speed (CAGESP), ball spin ( BSPEED ) , and ball spin speed 
(BALLSP). If other parameters as coolant inlet temperature, 
cage load, pressure drop, and coolant flowrate need to be 
varied, reference (14) states the affected inputs to 
SHABERTH/SINDA that must also be varied. These parameters 
will be considered fixed in this study. 

The SHABERTH inputs indicate a four-bearing system being 
modeled. However, due to the arbitrarily choosen small 
initial contact angle ©< 0 to be + 5 and zero diametrical 
clearance, bearings 1 & 2 are Hummy bearings in this model. 
Since the BSMT has four 57 mm bearings shown schematically 
in Figure 2, symmetry was used and only 2 of the 4 bearings 
are actually analyzed by SHABERTH. Therefore, bearings 1 
and 2 (the pump end bearings for HPOTP) are the dummy 
bearings and bearings 3 and 4 (the turbine end bearings for 
HPOTP) are analyzed. The SINDA model was written only for 
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bearing #3. The grid generation and nodal numbering was 
performed similiar to the process shown in (8,9 and 13) for 
the 45 mm pump-end bearings. The user need only be 
concerned with SINDA's coolant inlet and saturation 
temperatures (lines 697-709), cage heat (line 757), half of 
the viscous heat generation rates for bearing #4 (lines 760- 
761 for nodes 2 and 3) and for bearing #3 (lines 763-764 for 
nodes 5 & 6) and coolant flowrate per ball (lines 2228- 
2236). Also, specific heat vs. temperature lines 2293-2300 
of SINDA, must be changed when using different coolants. 
Notice in the initial nodal temperature guess in the 
SHABERTH input, only the 3rd line representing bearing #3 
has been deviated from an initial value of -170°F. These 
temperatures represent the shaft, inner ring, inner race, 
ball, outer race, outer ring, housing and fluid bulk 
temperatures. These initial temperatures will change with 
each iteration of SHABERTH/SINDA until either convergence or 
divergence occurs. Also, change the modulus of elasticity 
and thermal expansion coefficients to match the initial 
temperatures of bearing #3. They will also be updated in 
the iteration process. 

Axial preload can be included by setting the diametrical 
clearance of bearing #3 and #4 to a non-zero value. In the 
Appendices, a table is presented relating the amount of 
axial preload to the diametrical clearance. This was 
generated by running SHABERTH only at steady state 
temperature and denoting the Fx (x force reaction) in the 
output. Therefore, the amount of diametrical clearance 
inputted is related to the Fx force reaction which is the 
axial preload on bearings #3 and #4. These results are 
independent of coolant used and flowrate based on the 
simulation . 

In this study, the coefficient of friction was set at 0.2, 
tester pressure was 480 psia, saturated temperature for LN2 
was -233. 8°F and for LOX was -200. 8°F and the 
coolant flowrate was 6.4 lbm/sec. The axial preload -was set 
at 1000 lbs by setting the diametrical clearance input to 
be 0.013 mm on bearings #3 and #4. 


RESULTS 


Due to input parameter problems and UNIVAC down-time, 
the study of the BSMT axial load variation producing heat 
generation rates and contact stresses that were compared to 
a HPOTP simulation was abandoned at a shaft speed of 20,000 
rpm. At 30,000 rpm, a coolant flowrate of 4.6 lbm/sec was 
used initially for both the BSMT LN2 and LOX simulations. 
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This was the coolant flowrate used in the HPOTP simulation. 
At this flowrate for both LN2 and LOX coolants, the 
solutions diverged. The range of applied axial loads were 
from 1000 lbs to 3000 lbs with a fixed preload of 1000 lbs 
for these cases. As the axial load Q inc reased , the ball 
temperature accelerated toward 1000 F in 3 to 4 
iterations before divergence was declared. Based on these 
initial results, it was decided to increase the coolant 
flowrate to 6 . 4 lbm/sec for both LN2 and LOX BSMT models. 

In this process, however, several errors were found in the 
SINDA source data. Specifically, lines 2293-2300 were not 
changed in the LOX SINDA file. These lines list the 
specific h\eat vs. temperature of the coolant used. So, the 
LOX SINDA file was still using LN2 data. Also, in the 
SHABERTH input file, the LN2 fluid properties of specific 
heat, thermal conductivity, and Prandtl number had to be 
adjusted at the saturated temperatures. Since the tester 
operating pressure of 480 psi is near the critical pressure 
of LN2 of 493 psi, the variation in these properties were 
held at a constant value at the saturation temperature. 

This should stabilize the heat transfer conductance 
calculations according to SRS. So, these two problems could 
have played a part in the divergence of the solution at a 
flowrate of 4.6 lbm/sec. 

The above changes were made to the SHABERTH/SINDA input 
files and with the coolant flowrate value changed to 6.4 
lbm/sec, another series of program executions were 
performed. From this series of computer runs, Tables 1 and 
2 show the converged results of the heat generation rates 
and Hertz stresses in bearings 3 & 4. As shown, for both 
heat generation rates and Hertz stresses, there is no 
significant difference between using LN2 or LOX coolants 
for the range of axial loads. From Table 1, for bearing 
#3, there is a reasonable agreement between the BSMT and 
HPOTP simulations. For bearing #4, the BSMT simulation 
under predicts the HPOTP simulation by a factor of 1/2. 

This effect may be caused by the HPOTP simulation having a 
SINDA model of both bearings 3 & 4 whereas, the BSMT model 
only has bearing #3 thermally modelled. In Table 2, again, 
Inhere is no significant difference in Hertz stress for 
bearings #3 and #4 due to the coolant used in the BSMT 
model. The results from Table 2 show a reasonable 
agreement of outer and inner race Hertz stresses for the 
BSMT and HPOTP simulation for bearing #3; however, the BSMT 
model again underestimates Hertz stresses by about one- 
fourth compared to the HPOTP simulation. From these 
results, it is difficult to predict how much axial load 
could exactly predict the HPOTP simulation results. 

Further studies are necessary to attempt to establish an 
equivalent load relationship. 
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Table 1: Total Heat Generation Rates of the BSMT & HPOTP Simulations (Preload 1000 lbs) 

in Watts 
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Table 2: Hertz Contact Stresses of the BSMT & HPOTP Simulation (Preload-1000 lbs) 

in N/mm^ 



XXIII-10 


CONCLUSIONS AND RECOMMENDATIONS 


Based on my limited results, no relationship can be 
established at this time between the BSMT simulation and the 
HPOTP simulation loadings. In the BSMT simulation, no axial 
load above 1000 lbs (4446. 5N) would result in a stable 
thermally converging solution at a shaft speed of 30,000 
rpm. Based on this study, several recommendations for 
future research in this area are as follows. 

1. The continuation of this study at a lower shaft 
speed to determine it's effect on the comparison 
of heat generation rates contact stresses and 

on enabling the use of higher axial loads. 

2. The study of the effects of coolant flowrates and 
coefficient of friction on the comparison between 
BSMT simulation axial loads and HPOTP simulation 
loads . 

3. The investigation of other bearing parameters that 
need be included besides heat generation rates and 
contact stress in the equivalency of BSMT and HPOTP 
loading . 

4. The correlation of BSMT simulation axial load re- 
sults to actual BSMT tester data for both LN2 and 
coolants . 

Hopefully, from these recommendations, an equivalency 
between BSMT axial loads and HPOTP loadings can be found. 
However, the possibility exists that an applied axial load 
only may never produce equivalent HPOTP conditions in the 
bearing tester. So, the logical alternative may be to 
incorporate a workable radial load capability to the bearing 
tester and to the SHABERTH BSMT simulation. The alternative 
would lead to a matching of both axial and radial load 
conditions between the tester and turbopump to hopefully 
generate the same mechanical and thermal environment for 
the bearings. 

For SHABERTH's results to be a reliable predictor of 
bearing performance, it must have reliable inputs based upon 
both experimental data and analytical formulation. SHABERTH 
is also constantly being modified and updated by SRS to make 
it more versatile in its simulation of a shaft bearing system 
by including more bearing theory. Eventally, SHABERTH could 
become an important analytic tool for both the current HPOTP 
or BSMT configuration and for any future alternative 
configurations that may be developed. 
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APPENDICES 


SHABERTH/SINDA INPUTS 
FOR LN„ 
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lines 2296-2301 

specific heat vs temperature at 480 psia for LOX at node 27 


Diametrical Clearance vs Axial Preload 


Diametrical 


Axial 

Clearance 


preload 

Input to 


on 

SHABERTH 


bearing pair 

(mm) 

(N) 

( lb) 


0.0043 

0.009 

0.013 

0.0148 

0.025 

0.05 


5137 

4750 

4450 

4315 

3651 

2466 


1155.3 

1068.3 

1000.8 

970.4 

821.1 

554.6 
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